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Localization-induced coherent backscattering effect in wave dynamics 
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We investigate the statistics of single-mode delay times of waves reflected from a disordered wave- 
guide in the presence of wave localization. The distribution of delay times is qualitatively different 
from the distribution in the diffusive regime, and sensitive to coherent backscattering: The probabil- 
ity to find small delay times is enhanced by a factor close to \/2 for reflection angles near the angle 
of incidence. This dynamic effect of coherent backscattering disappears in the diffusive regime. 
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I. INTRODUCTION 

The two most prominent interference effects arising 
from multiple scattering are coherent backscattering and 
wave localization [|l|-||. Both effects are related to the 
static intensity of a wave reflected or transmitted by 
a medium with randomly located scatterers. Coherent 
backscattering is the enhancement of the reflected inten- 
sity in a narrow cone around the angle of incidence, and 
is a result of the systematic constructive interference in 
the presence of time- reversal symmetry Mbl. Localiza- 
tion arises from systematic destructive interference and 
suppresses the transmitted intensity Q. 

This paper presents a detailed theory of a recently 
discovered [|7| interplay between coherent backscattering 
and localization in a dynamic scattering property, the 
single-mode delay time of a wave reflected by a disordered 
waveguide. The single-mode delay time is the derivative 
(j)' — d(j)/Auj of the phase </> of the wave amplitude with 
respect to the frequency w. It is linearly related to the 
Wigner-Smith delay times of scattering theory [p|-p^ and 
is the key observable of recent experiments on multiple 
scattering of microwaves |0] and light waves |12]. Van 
Tiggelen, Sebbah, Stoytchev, and Genack ||l^ have de- 
veloped a statistical theory for the distribution of 0' in a 
waveguide geometry (where angles of incidence are dis- 
cretized as modes). Although the theory was worked out 
mainly for the case of transmission, the implications for 
reflection are that the distribution P(0') does not depend 
on whether the detected mode n is the same as the inci- 
dent mode m or not. Hence it appears that no coherent 
backscattering effect exists for P{(j)'). 

What we will demonstrate here is that this is true only 
if wave localization may be disregarded. The previous 
studies |ll|,|l^ dealt with the diffusive regime of wave- 
guide lengths L below the localization length ^. (The 
localization length in a waveguide geometry is ^ ~ Nl, 
with N the number of propagating modes and I the mean 
free path.) Here we consider the localized regime L > ^ 
(assuming that also the absorption length fa > C)- The 
distribution of reflected intensity is insensitive to the 
presence or absence of localization, being given in both 



regimes by Rayleigh's law. In contrast, we find that the 
delay-time distribution changes markedly as one enters 
the localized regime, decaying more slowly for large \4>'\- 
Moreover, a coherent backscattering effect appears: For 
L > £_ the peak of P{(j>') is higher for n ^ m than for 
n 7^ TO by a factor which is close to v2, the precise factor 
being V2 x ^§^ = 1.35. 

We also consider what happens if time- reversal symme- 
try is broken, by some magneto-optical effect. The co- 
herent backscattering effect disappears. However, even 
for n ^ m, the delay-time distribution for preserved 
time-reversal symmetry is different than for broken time- 
reversal symmetry. This difference is again only present 
for i > f and vanishes in the diffusive regime. 

The plan of this paper is as follows: In Section || we 
specify the notion [^ of the single- mode delay time </>', 
relate it to the Wigner-Smith delay times and review the 
results |13J for the diffusive regime, extending them to 
include ballistic corrections. This Section also contains 
the random-matrix formulation for the localized regime, 
that provides the basis for our calculations, and includes 
a brief discussion of the conventional coherent backscat- 



tering effect in the static intensity /. Section [II presents 
the calculation of the joint distribution of (p' and /. We 
compare our analytical theory with numerical simula- 
tions and give a qualitative argument for the dynamic 
coherent backscattering effect. The role of absorption is 
discussed, as well as the effect of broken time-reversal 
symmetry. Details of the calculation are delegated to the 
appendices. 



II. DELAY TIMES 

A. Single-mode delay times 

We consider a disordered medium (mean free path I) 
in a waveguide geometry (length L), as depicted in Fig. 
n. There are N ^ 1 propagating modes at frequency lu, 
given by A'^ = ttA/X^ for a waveguide with an opening of 
area A. The wave velocity is c, and we consider a scalar 
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mainly focus on the optically more relevant case of pre- 
served time-reversal symmetry, we will also discuss the 
case of broken time-reversal symmetry for comparison. 
These two cases are indicated by the index f3 = 1, 2, 
respectively. 



B. Relation to Wigner-Smith delay times 



FIG. 1. Sketch of a waveguide containing a randomly scat- 
tering medium and illuminated by a monochromatic plane 
wave. We study the frequency dependence of the phase of 
the reflected wave amplitude in a single speckle, correspond- 
ing to a single waveguide mode. The derivative cp' = d^/dw 
is the single-mode delay time. 



wave (disregarding polarization) for simplicity. In the nu- 
merical simulations we will work with a two-dimensional 
waveguide of width W, where N — 2W/X. 

We study the dependence of the reflected wave ampli- 
tude 



= Vle"* 



(1) 



on the frequency to. The indices n and m specify the 
detected and incident mode, respectively. (We assume 
single-mode excitation and detection.) Here / = |r„mP 
is the intensity of the reflected wave in the detected mode 
for unit incident intensity, and characterizes the static 
properties of the reflected wave. Dynamic information is 
contained in the phase derivative 



duj ' 



(2) 



which has the dimension of a time and is called the single- 
mode delay time [^l|^. The intensity / and the delay 
time (/)' can be recovered from the product of reflection 
matrix elements 



P = r„m{^ + \5uj)r*^„^{uj ~ \5uj) 



(3) 



evaluated at two nearby frequencies ijj±\5u). To leading 
order in the frequency difference 5uj one has 

Ini Q 

p^ I[l + i5uj(j)') ^ I = lim Rep, (j)' = lim -— ^ . 

(4) 

We seek the joint distribution function P(/, 0') in an 
ensemble of different realizations of disorder. We dis- 
tinguish between the diffusive regime where L is small 
compared to the localization length ^ ~ Nl, and the lo- 
calized regime where L ^ ^. Localization also requires 
that the absorption length ^a ^ C- We will contrast the 
case of excitation and detection in two distinct modes 
n ^ m with the equal- mode case n = ni. Although we 



In the localized regime (^ <C L,£^s) we can relate the 
single-mode delay time (j>' to the Wigner-Smith p-|ic|] 
delay times r^, with i = I, . . . , TV. The r^'s are defined for 
a unitary reflection matrix r (composed of the elements 
Tnm), hence they require the absence of transmission and 
of absorption. One then has 

= C/'i'diag(Ti,...,TAr)[/, (5a) 



dw 

dr^ 
—ir* — — = V'' diag (ti , 
duj 



,tn)V, 



(5b) 



with U and V unitary matrices of eigenvectors. In the 
presence of time-reversal symmetry r is a symmetric ma- 
trix, hence V = U in this case. 
For small Suj we can expand 

r{LU ± Ulu) = V'^U ± U Suj V'^ diag (n, . . . , rjv) [/ . (6) 



Inserting it into Eq. 
yields the relations 

0' = Re -i , / : 

^0 



and comparison with Eq. 



UnP 



A, 



E 



T, UiVi 



(7) 



We have abbreviated Ui = Uim, ^i — Vin- In the special 
case n = m the coefficients Ui and Vi are identical in the 
presence of time-reversal symmetry. 

The distribution of the Wigner-Smith delay times in 
the localized regime was determined recently p4| . In 
terms of the rates pi = l/r^ it has the form of the La- 
guerre ensemble of random-matrix theory. 



PiifJ'i}) oc n 1/^* ~ ^J 



Yl^it^k)e 



-7(,3Ar+2-,3)Mfc 



(8) 



i<j 



where the step function Q{x) — 1 for a; > and for 
X < 0. The parameter 7 is defined by 



7 = al/c, 



(9) 



with the coefficient a — 7r^/4 or 8/3 for two- or three- 
dimensional scattering, respectively. Eq. (j^ extends the 
iV = 1 result of Refs. (f|-0] to any iV. 

The matrices U and V in Eq. (||) are uniformly dis- 
tributed in the unitary group. They are independent for 
/3 = 2, while U = V ioi P = I. In the large- A^ limit the 
matrix elements become independent Gaussian random 
numbers with vanishing mean and variance 1/N. Hence 



= 0, 



\v,f)^N-\ 



(10) 



with Ui — Vi for n — m and (3 — 1. Corrections to this 
Gaussian approximation are of order 1/N. 



C. Diffusion theory 

The joint probability distribution P{I, (f)') in the diffu- 
sive regime I <^ L <^ S, has been derived in Refs. 0,0, 






/2)-l/2 



exp 



(11) 



Diffusion theory predicts that the distribution of delay 
times ( |l4|) as well as the values of the constants 0' and 
Q do not depend on the choice n — m ov n ^ m (and 
also not on whether time-reversal symmetry is preserved 
or not). Hence there is no dynamic effect of coherent 
backscattering in the diffusive regime. 



D. Ballistic corrections 



with constants /, 0', and Q. ft has the same form for 
transmission and reflection, the only difference being the 
dependence of the constants on the system parameters. 
Here we focus on the case of reflection, because we are 
concerned with coherent backscattering. 

From the joint distribution function ( |lf|) one obtains 
for the intensity the Rayleigh distribution 



1 



Pdiff(/) = jexp(-///). 



(12) 



Hence / is the mean detected intensity per mode. It is 
given by 



/ 



1 

TV 



(1 + 5[3i5n7n)^ 



(13) 



assuming unit incident intensity. The factor-of-two en- 
hancement in the case n = m is the static coher- 
ent backscattering effect mentioned in the introduction, 
which exists only in the presence of time-reversal sym- 
metry (/3 = 1). Eqs. dll) and ( p^ ) remain valid in the 
localized regime, since they are determined by scattering 
on the scale of the mean free path. Hence L ^ I \s suf- 
ficient for static coherent backscattering, and it does not 
matter whether L is small or large compared to ^. 

By integrating over / in Eq. (O) one arrives at the 
distribution of single- mode delay times [|rf],|3| , 

PM<t^')--^[Q+{^'/^'~m-'"'- (14) 

Hence (j)' is the mean delay time while \/Q sets the rel- 
ative width of the distribution. These constants are de- 
termined by the correlator pT|jr^ 



■ 7/ r 17/2/^ 



Diffusion theory gives 

4)' = 27s/3, Q = 2s/5. 
Here 7 is given by Eq. (||). We have defined 
s — a' L/l, 



(15) 
(16) 

(17) 



where the numerical coefficient 



2/7r, 3/4 for two- 



three-dimensional scattering. (The corresponding result 
for Q given in Ref. p3[ is incorrect.) 



The expressions for the constants /, 0', and Q given 
above are valid up to corrections of order 1/ L. Here we 
give more accurate formulas that account for these ballis- 
tic corrections. (We need these to compare with numer- 
ical simulations.) We determine the ballistic corrections 
for Q and 0' by relating the dynamic problem to a static 
problem with absorption. (This relationship only works 
for the mean. It cannot be used to obtain the distribution 
|l9|.) The mean total reflectivity 



a — 1 + X — V 2a; + x^ coth s\/2x + x'^ + arcosh (I + x) 

(18) 

for absorption a'x per mean free path has been evaluated 
in Ref. [gOJ. [Here a' is the same constant as in the 
definition of s, Eq. (|l|).] We identify C12 = a{x)/a{0) 
by analytic continuation to an imaginary absorption rate 
X = —iSujj. Expanding in x to second order we find 



0' =7 



s{3 + 2s) 
3(1 + s) ' 



Q 



8s3 + 28s2 + 30s + 15 
5(2s + 3)2 ■ 



(19) 



E. Numerical simulation 

The validity of diffusion theory was tested in Refs. 
||ll|-[l3[ by comparison with experiments in transmission. 
In Fig. we show an alternative test in reflection, by 
comparison with a numerical simulation of scattering of a 
scalar wave by a two-dimensional random medium. (We 
assume time- reversal symmetry.) The reflection matri- 
ces r{uj ± ^Suj) are computed by applying the method 
of recursive Green functions |21|] to the Helmholtz equa- 
tion on a square lattice (lattice constant a). The width 
W = 100 a and the frequency uj — 1.4 c/a are chosen 
such that there are TV = 50 propagating modes. The 
mean free path I — 14.0 a is found from the formula |22] 
trrr^ = Ns{l + s)^^ for the reflection probability. The 
corresponding localization length <^ = NL/s — 1100 a. 
The parameter 7 = 46.3 a/c is found from Eq. ( |l9| ) 
by equating {<j)') = 0'. (This value of 7 is somewhat 
bigger than the value 7 = -k^I/Ac = 34.5 a/c expected 
for two-dimensional scattering, as a consequence of the 
anisotropic dispersion relation on a square lattice.) We 
will use the same set of parameters in the interpretation 
of the results in the localized regime, later in this paper. 




FIG. 2. Distribution of the single-mode delay time 4>' in 
the diffusive regime. The result of numerical simulation (data 
points) with A'' = 50 propagating modes is compared to the 
prediction ( |l4| ) of diffusion theory (solid curve). There is no 
difference between the case n = m oi equal-mode excitation 
and detection (open circles) and the case n 7^ m, of excitation 
and detection in distinct modes (full circles). 



Our numerical results confirm that in the diffusive 
regime the distribution of delay times (/)' does not distin- 
guish between excitation and detection in distinct modes 
(n 7^ m, full circles) and identical modes {n = m, open 
circles) . 



III. DYNAMIC COHERENT BACKSCATTERING 
EFFECT 

A. Distinct-mode excitation and detection 

We now calculate the joint probability distribution 
function P{I, (f)') of intensity / and single-mode delay 
time (f)' in the localized regime, for the typical case n ^ m 
of excitation and detection in two distinct modes. We as- 
sume preserved time-reversal symmetry (/? = 1), leaving 
the case of broken time-reversal symmetry for the end of 
this section. 

It is convenient to work momentarily with the weighted 
delay time W — (j>'I and to recover P{I, 4>') from P{I, W) 
at the end. The characteristic function 



x{p, q) = { 



^ — ipI — iqW 



) 



(20) 



is the Fourier transform of P{I, W). The average (• • •) is 
over the vectors u and v and over the set of eigenvalues 
{ti}. The average over one of the vectors, say v, is easily 
carried out, because it is a Gaussian integration. The 
result is a determinant, 



X{p,q) ^{det{l + iH/N)-^) , 

TT * T I 1 /-* T I *-T\ 

H ~ pu u + ^^(u u 4- u u ) 



(21a) 
(21b) 



The Hermitian matrix iJ is a sum of dyadic products 
of the vectors u and u, with Ui = UiTi, and hence has 



only two non- vanishing eigenvalues A_|_ and A_ . Some 
straightforward linear algebra gives 



A± = 5 [qBi +p± y/2pqBi + q^B2 + p^ 
where we have defined the spectral moments 



b.-y: 



u,\\f 



(22) 
(23) 



The resulting determinant is 

det(l + H/N)-^ = (1 + \+/N)-^{l + \-/N)-^, (24) 
hence 

„2 



x(p, q) = 



1 



■ip iq 

N N 



Bi 



47V2 



(B2 - Bi) 



(25) 



An inverse Fourier transform, followed by a change of 
variables from /, W to /, 0', gives 

p{i,(t>') = e(/)(iv3j/^)i/2g-A'/ 



(B2-5?)"^/^exp( -iV/ 



B2 - Bl 



(26) 



The average is over the spectral moments Bi and i?2, 
which depend on the UiS and r^'s via Eq. (23). 

The calculation of the joint distribution P{Bi,B2) is 
presented in Appendix ^. The result is 



P(Bi,S2)-e(Si)e(B2)exp - 



Bf-fN^ 
^2 



{B2+-fN^Bi)e^p 



NBl 
2-/N\ 



^-^(2^2 

4B| ^ 2 



ABIB2N + BfN^)Ei 



2jN 

"bT 



(27) 



where Ei (x) is the exponential-integral function. The 
distribution P(/, (/)') follows from Eq. ( p6| ) by integrating 
over Bl and B2 with weight given by Eq. (27). 



Irrespective of the distribution of Bi and B2, we re- 
cover from Eq. (p3) the Rayleigh law (W% for the inten- 
sity /. The distribution P{(j)') = J^ dl P{I,(j)') of the 
single-mode delay time takes the form 



P(0') = / dPi / dB 



P{Bi,B2){B2^Bl) 
2{B2 + <j)'^ - 2Bi(j)'f/'^ 



(28) 



In Fig. y this distribution is compared with the result of a 
numerical simulation of a random medium as in Section 
|IIE| , but now in the localized regime. The same value 
for 7 was used as in Fig. 0, making this comparison a 
parameter-free test of the theory. (Note that 7 alone 
determines the complete distribution function in the lo- 
calized regime, in contrast to the diffusive case where 
two parameters are required.) The numerical data agrees 
very well with the analytical prediction. 
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FIG. 3. Distribution of the single-mode delay time cj)' in 
the localized regime. The results of numerical simulations 
with A'^ = 50 propagating modes (open circles for n — m, full 
circles for n ^ m) are compared to the analytical predictions. 
The curve for different incident and detected modes n ^ m 
is obtained from Eqs. (|27[ ) and (|2q ). The curve for n = m, is 
calculated from Eqs. (|2£|) and (|30p. The same value for 7 is 
used as in the diffusive regime, Fig. 0. 
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FIG. 4. Distribution of the single-mode delay time </!>' in 
the localized regime for preserved time-reversal symmetry, in 
the limit A'' — + oq. In this limit P{4>') becomes symmetric for 
positive and negative valu es of 0'. Compared are the result 
for n 7^ m [Eqs. (H), ( |A16| )] and n = m [Eqs. g, @, (||)]. 
The distribution for n = m falls on top of the distribution for 
n ^ m when (f)' is rescaled by a factor 1.35 (dashed curve, 
almost indistinguishable from the solid curve for n 7^ m). 



B. Equal-mode excitation and detection 

We now turn to the case n = to of equal-mode ex- 
citation and detection, still assuming that time-reversal 
symmetry is preserved. Since Ui = Vi, we now have 



Re 



Co 



\Co\ 



Ck 



E^'"'- (29) 



The joint distribution function P(Co,Ci) of these com- 
plex numbers can be calculated in the same way as 
P{Bi, ^2). In Appendix |c| we obtain 

P(Co,Ci)cxexp(-A^|Co|V2) / dss^e'' 

Jo 

x(l + ^-|^ReC„C,y\ (30) 

The corresponding distribution function P{(j)') is plotted 
also in Fig. and compared with the results of the nu- 
merical simulation. Good agreement is obtained, without 
any free parameter. 



C. Comparison of both situations 

Comparing the two curves in Fig. 0, we find a striking 
difference between distinct-mode and equal-mode excita- 
tion and detection: The distribution for n — m displays 
an enhanced probability of small delay times. In the 
vicinity of the peak, both distributions become very sim- 
ilar when the delay times for n ^ m are divided by a 
scale factor of about \/2. In limit N -^ 00 (see following 



sub-section), the maximal value Pi(t>Lca.k) ~ \/2/nN^j'^ 
for n = TO, is larger than the maximum of P{(j)') for n ^ m 
by a factor 



P. 



peak/ 



P, 



= V2: 



n^m.y'Pponk) 



4096 
137l7r 



= 1.35. 



(31) 



Correspondingly, the probability to find very large delay 
times is reduced for n = m. This is reflected by the 
asymptotic behavior 



P(0') 



7iV3/2 



(27r)-i/2 
y^/4 



for n = m, 
for n ^ m. 



(32) 



The enhanced probability of small delay times for 
n = TO is the dynamic coherent backscattering effect 
mentioned in the introduction. The effect requires lo- 
calization and is not observed in the diffusive regime. 



D. Limit iV ^ 00 

The results presented so far assume A^ ;:^ 1, but retain 
finite- iV corrections of order N"^'^. (Only terms of order 
1/A^ and higher are neglected.) It turns out that the 
asymmetry of P{4>') for positive and negative values of 
(p' is an effect of order N~^/'^. The asymmetry is hence 
captured faithfully by our calculation. We now consider 
how the asymmetry eventually disappears in the limit 
iV^ 00. 

For distinct modes n ^ ni, the spectral moments scale 
as Bi ~ 7A^ and B2 ~ 7^A^^. With <f)' ~ jN^/^, one 
finds that Bi can be omitted to order N^^^"^ in Eq. (|2q). 
One obtains the symmetric distribution 



r 

Jo 



dBo 



P{B2)B2 



(33) 



2(^2 + 0'2)3/2 ' 

plotted in Fig. |. 

For identical modes n = m, observe that the quantities 
Co and Ci become mutually independent in the large- A^ 
limit: The cross-term (77V) "^ Re Co C^* in Eq. @ is of 
relative order N^^/"^ because Co ~ N^^^^ and Ci ^ 
^N. Hence to order N^^^^ the distribution factorizes, 
P(Co,Ci) = P(Co)P(Ci). The distribution of Co is a 
Gaussian, 



P(Co) = ^exp(-iV|Co|V2), 

as a consequence of the central-limit theorem, and 

-5/2 



P(Ci)cx / dss^e^' [1 + 



|CiPs2 



(34) 



(35) 



^2^2 

The resulting distribution of cj)' = Re (Ci/Co) is also plot- 
ted in Fig. |. 

The dynamic coherent backscattering effect persists in 
the limit A^ ^- 00, it is therefore not due to finite- A^ 
corrections. The peak heights differ by the factor given 
in Eq. (II. 

E. Interpretation in terms of large fluctuations 

In order to explain the coherent backscattering en- 
hancement of the peak of P{(f)') in more qualitative terms, 
we compare Eq. (E9h for n = m with the corresponding 
relation (|^) for n ^ m. 

The factorization of the joint distribution function 
P(Co,Ci) discussed in the previous sub-section can be 
seen as a consequence of the high density of anoma- 
lously large Wigner-Smith delay times r^ in the La- 
guerre ensemble (&). The distribution of the largest time 
Tmax — maxi Ti follows from the distribution of the small- 
est eigenvalue in the Laguerre ensemble, calculated by 
Edelman |23]. It is given by 



^(Tmax) = -^ exp(-7iV^/r,nax)- 



(36) 



As a consequence, the spectral moment Ci is dominated 
by a small number of contributions uf r^ (often enough 
by a single one, say with index i = 1), while Cq can 
be safely approximated by the sum over all remaining 
indices i (say, i ^ 1). The same argument applies also 
to the spectral moments Ak which determine the delay- 
time statistics for n ^ m, hence the distribution function 
P{Aq,Ai) factorizes as well. 

The quantities Aq and Co have a Gaussian distribution 
for large A^, because of the central-limit theorem, with 
P(Co) given by Eq. <M) and 



P(Ao) = -exp(-A^|Aon. 

TT 



(37) 







FIG. 5. Distributions of Bi = Bif-yN and 62 = Bi/'y^N^. 
The a naly tic pre dicti on from Eq. ( P7|) [for explicit formulas see 
Eqs. (A15) and (A16)] is compared to the result of a numerical 
simulation of the Laguerre ensemble with TV — 50. 



It becomes then clear that the main contribution to the 
enhancement ( |3l| ) of the peak height, namely the factor 
of v^, has the same origin as the factor-of-two enhance- 
ment of the mean intensity /. More precisely, the relation 
P{Ao ^ x) = 2 P(Co ^V2x) leads to a rescaling of P(/) 
for n = 771 by a factor of 1/2 and to a rescaling of P{4)') 
by a factor of -\/2- The remaining factor of ^^^i^ = 0.95 
comes from the difference in the distributions P{Ai) and 
P(Ci). It turns out that the distribution 



47r7Ar (4-h \Ai/-iN\^s^f 



P{Ai) - / As 
Jo 

X [e""(64 + 32s + Us^ + s^) - 3s^Ei (s)] (38) 

(derived in Appendix |D|) is very similar to P[Ci) given 
in Eq. (pq), hence the remaining factor is close to unity. 
The large r^'s are related to the penetration of the 
wave deep into the localized regions and are eliminated 
in the diffusive regime i ^ ^. In the following sub-section 
we compare the localized and diffusive regimes in more 
detail. 

F. Localized versus diff'usive regime 



Comparison of Eqs. (|ll|) and (|26| ) shows that the two 
joint distributions of / and </>' would be identical if statis- 
tical fluctuations in the spectral moments Pi, P2 could 
be ignored. The correspondences are 



Pi 



P2 — Pi 



(39) 



However, the distribution P(Pi,P2) is very broad (see 
Fig. ph, so that fluctuations can not be ignored. The 
most probable values are 



P 



typical 



7iV, B. 



typical 



YN^ 



(40) 



but the mean values (Bi), {B2) diverge — demonstrating 
the presence of large fluctuations. In the diffusive regime 
L < ^ the spectral moments Bi and B2 can be re- 
placed by their ensemble averages, and the diffusion the- 
ory |0,f3[ is recovered. (The same applies if the absorp- 
tion length ^a ^ ^0 

The large fluctuations in Bi and B2 directly affect the 
statistical properties of the delay time 0'. We compare 
the distribution (p8) in the localized regime (Fig. 0) with 
the result (|lj) of diffusion theory (Fig. 0). In the lo- 
calized regime the value 0' ^^ — B^'^^'^'^ at the center 
of the peak of P{4>') is much smaller than the width of 
the peak A0' ~ (B^^P'^^y/^ ~ (ji'^.^^i^/lf^. This holds 
also in the diffusive regime, where </>' ^k = 4>' ^-nd A0' ~ 
'^poak(^/0^''^- However, the mean {4>') = {Bi) diverges 
for P, but is finite (equal to 0') for P^m- For large B2 

-rNr'I'^^B:^^!'^. As a 



one has asymptotically P(B2) 
consequence, in the tails P{(j>') falls off only quadratically 
[see Eq. (p2h], while in the diffusive regime Pdis{4'') ^ 
^Qcj)' 10'!^^ falls off with an inverse third power. 



G. Role of absorption 

Although absorption causes the same exponential de- 
cay of the transmitted intensity as localization, this de- 
cay is of a quite different, namely incoherent nature. The 
strong fluctuations in the localized regime disappear as 
soon as the absorption length ^a drops below the local- 
ization length ^, because long paths which penetrate into 
the localized regions are suppressed by absorption. In 
this situation one should expect that the results of diffu- 
sion theory are again valid even for L }^ ^. This expecta- 
tion is confirmed by our numerical simulations. (We do 
not know how to incorporate absorption effects into our 
analytical theory.) 

In Fig. y we plot the delay-time distribution for two 
values of the absorption length Ca < C ^^'^ o^^ value 
Ca > Ci both for equal- mode and for distinct-mode ex- 
citation and detection. The length of the waveguide 
is L = 4.1^. The result for strong absorption with 
^a = 0.11^ is very similar to Fig. 0. Irrespective of 
the choice of the detection mode, the data can be fitted 
to the prediction (O) of diffusion theory. The plot for 
^a = 0.47i^ shows that the dynamic coherent backscatter- 
ing effect slowly sets in when the absorption length be- 
comes comparable to the localization length. The data 
also deviates from the prediction of diffusion theory. The 
full factor ( pi| ) between the peak heights quickly devel- 
ops as soon as ^a exceeds ^, as can be seen from the 
data for ^a = 2.1^. Moreover, this data can already be 
fitted to the predictions of random-matrix theory, with 
7 « 53.2 a/c. (The value 7 = 46.3 a/c of Sec. |IE| is 
reached when absorption is further reduced.) 
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FIG. 6. Single-mode delay-time distribution P{<l>') in the 
presence of absorption. The data points are the result of a nu- 
merical simulation of a waveguide with length L = 4.5^. Open 
circles are for equal-mode excitation and detection n = m, 
full circles for the case of distinct modes n ^ m. In the upper 
panel (with ^a < Oi the data is compared to the prediction 
( [L4|) of diffusion theory. In the lower panel we compare with 
the predictions ( pTJ - pOl ) of random-matrix theory. 
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FIG. 7. Comparison of the single-mode delay-time dis- 
tributions for preserved and broken time-reversal symmetry. 
The number of propagating modes is A*' = 50. The curves are 
calculated from Eq. (|8|) with P{Bi,B2) given by Eq. (127 
(/3 = l)orEq. (gT]) (/3 = 2). 



H. Broken time-reversal symmetry 

The case /3 = 2 of broken time-reversal symmetry is 
less important for optical applications, but has been re- 
alized in microwave experiments ||24|-t26|. There is now 
no difference between n = m and n ^ m. The matrices 
U and V have the same statistical distribution as for the 
case of preserved time-rever sal sy mmetry. Hence, by fol- 
lowing the steps of Section [II A , we arrive again at Eq. 
(p6|), with spectral moments B^ as defined in Eq. (23). 
Their joint distribution has now to be calculated from 
Eq. (pi) with [3 = 2. This calculation is carried out in 
Appendix ^. The result is 



P{Bi,B2) = ^^^3^' exp{~NBf/B2 - 2jN/Bi). (41) 

The distribution of single- mode delay times P{(f>') is given 
by Eq. (|8|) with this new function P{Bi,B2). We plot 
P{4>') in Fig. and compare it to the case of preserved 
time-reversal symmetry. The distribution is rescaled by 
about a factor of two towards larger delay times when 
time-reversal symmetry is broken. This can be under- 
stood from the fact that the relevant length scale, the lo- 
calization length, is twice as big for broken time-reversal 
symmetry (^ = 2NL/ s, while S, = NL/s for preserved 
time-reversal symmetry) . 



IV. CONCLUSION 



effect in the reflected intensity. The dynamic effect can 
be understood from the combination of the static effect 
and the large fluctuations in the localized regime. 

In the diffusive regime there is no dynamic coherent 
backscattering effect: The distribution of delay times is 
unaffected by the choice of the detection mode and the 
presence or absence of time-reversal symmetry. The ef- 
fect also disappears when the absorption length is smaller 
than the localization length. In both situations the large 
fluctuations characteristic for the localized regime are 
suppressed. 

Existing experiments on the delay-time distribution 
pd] , p^ have verified the diffusion theory jl^] . The theory 
for the localized regime presented here awaits experimen- 
tal verification. 
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APPENDIX A: JOINT DISTRIBUTION OF Bi 

AND Ba FOR /3 = 1 

We calculate the joint probability distribution function 
P (-61,-82) of the spectral moments Bi and B2, defined 
in Eq. ([2^), which determine P{I, (f)') from Eq. (|2^). We 
assume preserved time-reversal symmetry (/? = \)7 Since 
Bk = ^. JMip/i" , we have to average over the wavefunc- 
tion amplitudes w^, which are Gaussian complex numbers 
with zero mean and variance 1/N , and the rates /z; which 
are distributed according to the Laguerre ensemble (||) 
with (3=1. This Laguerre ensemble is represented as 
the eigenvalues of an iV x iV Hermitian matrix W^W, 
where 14^ is a complex symmetric matrix with the Gaus- 
sian distribution 



P{W) (X exp[--f{N + l)trW^W]. 



(Al) 



The calculation is performed neglecting corrections of or- 
der 1/iV, so that we are allowed to replace N + I hy N. 
The measure is 

dVF = Jl d Re Wij d Im Wij J| d Re Wu d Im Wu . ( A2) 

i<j i 



We have presented a detailed theory, supported by nu- 
merical simulations, of a recently discovered 0] coherent 
backscattering effect in the single-mode delay times of a 
wave refiected by a disordered waveguide. This dynamic 
effect is special because it requires localization for its ex- 
istence, in contrast to the static coherent backscattering 



1. Characteristic function 

In the first step we express P{Bi, B2) by its character- 
istic function. 



PiBi,B2) = 



(2^)2 



N 



dp 



(A3) 



x{p, 9) == ( n °^p 



\i=i 






(A4) 



and average over the m;'s, 



N 



\l=l 



x{p,q) = {l[(i + i^ + ^ '^ 



det (VKtVK)2 



\ det [(Tytw^)2 + ip{W^W)/N + iq/N] 

(A5) 

We have expressed the product over eigenvalues as a ra- 
tio of determinants. We write the determinant in the 
denominator as an integral over a complex vector z, 

x{p,q)oc I AW /dzexp[-7ArtrW^tvi/]det(VF^l^)2 

X exp{-z^[(VKtVF)^ + ip{W'^W)/N + iq/N]^} . (A6) 
This integral converges because W^W is positive definite. 

2. Parameterization of the matrix W 

Now we choose a parameterization of W which facili- 
tates a stepwise integration over its degrees of freedom. 
The distribution of W is invariant under transformations 
W -^ U^WU, with any unitary matrix U. Hence we 
can choose a basis in which z points into direction 1 and 
write W in block form 

^=(xa)- (A^) 

Here a is a complex number. For any N — 1 dimensional 
vector X we can use another unitary transformation on 
the X block after which x points into direction 2. Then 
W is of the form 

/ a X 0^ \ 
W^\ X b y^ \, (A8) 

\o y y J 

with the real number x — |x|. In this parameterization 

{W^W)n^\a\^ + x\ 



[{W^W)\ 

det W = [a{b - y^y^ V) - a^^ldet Y 



i\a\^+xY+x^j/ + x^\a + b*\^ 
et 
2y^ + tiY 



trW^'W = |ap + \b\^ + 2x^ 
dVF^d^ad^fedxdydr , 

with 2/ = |y|. A suitable transformation on Y allows to 
replace the term y^F~^y by y^{Y~^)ii. 

For this parameterization of W, the integrand in Eq. 



(A6) depends on the vectors x, y, and z only by their 
magnitudes x, y, and z = |z|. Hence we can replace 



dx 



p2Af-3 



da;, dy 



y 



2N-5 



dy, and dz 



z^^-^dz. 



3. Integration 



The integrand in Eq. (A6) involves z, p, and q in the 
form 



-z^{[{W''W)%i + ip{W^W)n/N + iq/N)] . (A9) 



cxp 



It is convenient to pass back to P{Bi,B2) by Eq 
because integration over p and q gives delta functions 



6[Bi - z2(|a|2 + x^)/N]6{B2 - z^N) 

= ■^S{B,/B2 - lap - x^)6iB2 ~ z^/N) . 

t>2 

Subsequent integration over z results in 

P(Bi,S2)cx I d^ad%dxdyx^^-^y^^-^ 

xB^-^S{Bi/B2-\a\^-x^) 
X [co\ab - a;2|4 + Ac2\a\''y^\ab - x^\^ + Ci\a\^y^] 
■Bl 



(AlO) 



X exp 



'NB2[ -^+x''y'' + x'\a* +b{- 



2-fNy^ 
(All) 



Here we omitted a term 77V(|ap -I- |6p -t- 2x^) in the ex- 
ponent, because it is of order 1/N as we shall see later. 
Furthermore, we denoted 



{\det Y\^\{Y-^)nn 
(|dety|4) 



(A12) 



These coefficients will be calculated later, with the re- 
sult Co = 1, C2 = 27, C4 = 47^. Integration over 
y yields for the terms proportional to Cm the factors 
{B2X^ + 27)""^"^+^, which can be combined with the 
factor {B2x'^)^~^, giving to order 1/A^ [we anticipate 
-f/B2X^ - Oil/N)] 



{B2X 



2\N-2 



{B2X^ + 27)^-2+™ 



{B2X 



2\-m 



exp 



27iV\ 
B^J 



(A13) 



We introduce a new integration variable by 6' = 6 + o* 
So far P{Bi, B2) is reduced to the form 



P{Bi,B2)(x I d^ad^b'dxx5{BilB2 

4c2|a|2 



lb' - 



X exp 



El 
B2 

2-fN 



lb' 



Bj_ 
B2 



\af-x') 



Ci\a\' 



Blx^ 



N^-NB2X^\b'\' 
B2 



B^x^ 



(A14) 



Let us now convince ourselves with this expression that 
we were just ified to omit the te rm "f N[\a\^ + \b\^ + 2x'^) 
in Eq. (All) and to use Eq. (A13). Indeed, the var- 
ious quantities scale as Bi ~ 7A, B2 — ^"^N^, and 



\b\' 



1/jN'^, because any 7 and N depen- 



dence disappears if one passes to appropriately rescaled 
quantities Bi/jN, etc. The terms omitted are therefore 
of order 1/A^. 



The remaining integrations in Eq. (A_14) are readily 
performed, with the final result (P7|). The distribution of 
Bi to order 1/N is 



PiBi) = ^{B, 



-Bf 



27iV)expf-^) 



(A15) 



The spectral moment Bi appeared before in a different 
physical context in Ref . [ po| , but only a heuristic approx- 
imation was given in that paper. Eq. (A15) solves this 
random-matrix problem precisely. 

For completeness we also give the distribution of the 
other spectral moment B2 (rescaled as B2 — -827^^^ '^) 
in terms of Meijer G functions, 



^="(^2) = 



1 



64S. 



5/2 



.1/2 



147r^2-16G[;';](B2| 



1 ^^ 

2' "' 2' 



20G°;;](B2|-iil) + 22G°;;](S2|i,i,l) 



^0,3/ R 



8G'3;o(-B2l5,l,|) +4:Gj'q{B2\^, |,2) 



'^4,1 



3G 



4,1 



B2 
B2 



—- n - 2 
1 

113 3 
2' 2' 2' 2 





-16G. 



1.3 

4,1 



B. 



n i - 2 

"' 2' 2'^ 
1 



(A16) 



4. Coefficients 



No w we calculate the coefficients C2 and C4 defined in 
Eq. (A12). It is convenient to resize the matrix Y to 
dimension TV (instead of iV — 2) and to set momentarily 
-fN — 1. We use again a block decomposition. 



Y 



a w 
w Z 



and employ the identities 



det Y = {a- w'^ Z~^w)det Z 
(r-i)n = (a - w^ Z-^w)-^ . 



Hence 



C4 = 



(|detZ|- 



(|detr|4) iN+l){N + 3) 
where we used Selberg's integral E^ for 

1 r(iv -F 4)r(iv + 2) 



(Idetrn^ 



6 



22JV 



In order to evaluate 

(|detZ|4(|ap 



C2 



/^Z-lw|2)) 



(|dety| 



(A17) 



(A18a) 
(A18b) 



(A19) 



(A20) 



(A21) 



it is again profitable to use unitary invariance and turn 
w into direction 1, 



w'Z-'w\'=w''\{Z-')ii\^ 



(A22) 



From (w^) = jN{N + 1) and (|ap) == 1 we obtain then 
the recursion relation 



C2{N) 



N 



iN + l){N + 3) N + 3 
which is solved by 

C2{N) 



N + 1 



C2{N - 1) , (A23) 



(A24) 



In order to reintroduce 7 we have to multiply c^ by 
(77V)™/2, and obtain to order 1/iV 



C2 = 27 , C4 = 47 

as advertised above. 



(A25) 



APPENDIX B: JOINT DISTRIBUTION OF Bi 

AND B2 FOR l3 = 2 

For broken time-reversal symmetry, the distribution 
of Bi and B2 has to be calculated from the Laguerre 
ensemble (||) with (3 — 2. Similarly as for preserved time- 
reversal symmetry, this ensemble can be obtained from 
the eigenvalues of a matrix W^W. The matrix W is once 
more complex, but no longer symmetric (it is also not 
Hermitian). It has the Gaussian distribution 



P{W) ex exp(-277VtrT4^^Ty), 
with measure 

dl^ = [| d Re Wij d Im Wij . 

id 



(BI) 



(B2) 



It is instructive to calculate first P{Bi), because it will 
be instrumental in the calculation of P{Bi,B2)- After 
averaging over the Ui's, the characteristic function takes 
the form 



X{p) = (exp(-ipSi)) = 



det W^W 



det {WW + ip/N) 



(B3) 



We express the determinant in the denominator as an 
integral over a complex vector z. Due to the invariance 
W -^ UWV of P{W) for arbitrary unitary matrices U 
and y, we can turn z into direction 1 and write 



W^ 




(B4) 



Then 
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X (|a|2 + d2X^x'^) exp[-(z2 + 2-fN){\a\'' + x^)] 
X exp[ip(Bi - z^/N) - 27iVa;'2] . (B5) 



Selberg's integral |Q gives 

(|detr|2|(y-i)np) 27iV 



rf2 = 



(|dety| 



7V-1 



(B6) 



The integration over p gives (5(z^ — NBi) and allows to 
eliminate z. The integration over x' amounts to replacing 



i;'2 = {N - 1)/277V 



The final integrations are 



most easily carried out by concatenating a to x, giving 
an A'^-dimensional vector y. Then 



P{Bi) « 



Ayy''''+^B'^-^eM~N{Bi 



27)r 



(xBf-XSi + 27)- 
which to order 1/A^ becomes 



-N-l 



(B7) 



P(i?i) = ^exp(-277V/i3i). 



The first steps in the calculation of the joint distribu- 
tion function of i?i and B2 are identical to what was done 
in Appendix H, and result in t he characteristic function 
x(p, 5) in the form of Eq. (A6), but with 7 replaced by 
27. Due to the unitary invariance of the W ensemble we 
can write 



/ a 



W 



\ 



X 

X b 
y 




0^ \ 
' 0^ 



y 



Y 



(B9) 



/ 



One integrates now over p and q and obtains delta func- 
tions as in Eq. ( AlO ). This is followed by integration over 
z. The calculation is then much simplified by recognizing 
that one can rescale the remaining integration variables 
in such a way (namely by introducing a^ = cP'Bi/B2, 



x^=x^B,/B2,y'^ 



~/2 - 

y X 



^B-^) that 



P{BuB2) = B^''cM-NBf/B2)f{Bi). 



(BIO) 



It is not necessary to give here f{Bi) as a lengthy multi- 
dimensional integral, since its functional form is easily 
recovered from the relation 

P(Si) = fdB2 P(Si, B2) = N-^B^^f{Bi). (Bll) 



We compare this with Eq. (|B8| ) and arrive at Eq. ( |4l[ ) 
The distribution of B2 has the closed-form expression 



P{B2) - ^-''N-^G"/o{j-^N-''B2 



1 

2' " 



(B12) 



APPENDIX C: JOINT DISTRIBUTION OF Co 

AND Ci 

We seek the joint distributions of the spectral moments 
Co and Ci, which determine (j)' and / for (3 — 1 and 
n = m via Eq. (E9). We start with the characteristic 
function 



x(Po,Pi) = (exp[iRe(poCo +P1C1)]) 



(CI) 



where po and pi are complex numbers, as are the quanti- 
ties Co and Ci themselves. Since Ck = J2i '^'i'''i ^^ have 
to average over the t^'s and the u^'s. Averaging over the 
u,'s first, we obtain 



xiPo,Pi) 



n 



1 



\pin +Po\ 

7V2 



-l/2\ 



(C2) 



We regard again the rates ^i = t~^ as the eigenvalues of 
a matrix product YY\ where Y will be specified below. 
Then the product of square roots can be written as a 
ratio of determinants, 



(B8) nfi 



\pin+po? 



7V2 



X det 



T\2 



(YY^) 



N' 



-1/2 



bol 



N^ 



det YY^ 






\Pi 



21-1/2 



7V2 



(C3) 



We will express the determinant in the denominator 
as a Gaussian integral over a real A'^-dimensional vector 
z. Hence it is convenient to choose Y real as well, so 
that one can use orthogonal invariance in order to turn 
z into direction 1. Moreover, there is a representation 
of Y which allows to incorporate the determinant in the 
numerator into the probability measure: We take y as a 
rectangular N x {N + 3) matrix with random Gaussian 
variables, distributed according to 

P{Y) ex exp(-77Vtrrr^) . (C4) 

The corresponding distribution of the eigenvalues /ii of 
YY^ is given in Ref. |g^ and differs from the Laguerre 
ensemble (|8|) by the additional factor Hi Mi = det YY^. 
In this representation, 

Xoc I dzz''-'(^exp{-z'{l + \p,\yN')[(YY^)']n} 
'2Repopl 



X exp 



iV2 



-[YY^U 



\Pi\' 



m 



(C5) 



where the average is now over Y . Inverse Fourier trans- 
fornration with respect to po and pi results in 



P(Co,Ci)(x 



,^-5 



eM-A{YY^f]ii] 



X exp 



^2K.-([rr^]n)2 

r2 \r< \wT^ n |2 



[(rr^)2jn 



\c^\'N' N' |Co-[yr^]iiCi| 



4z2 4^2[(yyT)2]^^_([yyT]^^)2 



(C6) 
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The orthogonal invariance of YY'^ allows us to param- 
eterize Y as 



Y 




(C7) 



with real numbers w > 0, w > 0, y > 0, a, and b, and an 
(A'' — 1) X {N + 1) dimensional matrix Z. It is good to 
see that Z drops out of the calculation, because it does 
not appear in 

[YY%i^a^ + v\ (C8a) 

[{YY^f]n = (a' + v^f + {aw + vbf + v^y' . (C8b) 

We replace b — b' — aw/v and introduce z' — zyv. The 
integral over z' can be written in the saddle-point form 
/dz'z'^e-^'V(^') oc f{y/N/2) for large N. The re- 
sulting expression varies with respect to the remaining 
variables on the scales 



N^a 



3„2 



N%'^ 



N^v^ ~ Ny^ 



w 



Oh-'). (C9) 



We use the given orders of magnitude to eliminate 
terms of order iV"^, but keep the residual correlations 
ReCaCl/-/N = 0{N-^^^). The joint distribution func- 
tion of Co and Ci is then 



P(Co,Ci) ex I dadb'dvv^dww^-^dyyexpi- 



lNy' 



X exp 
X exp 



-iNw'{i + ^^)~^{v^ + y' + b'') 



Nv^ReCoCl- 



Nv'y'\Ci\^ 



A^ICnP 



(CIO) 



2 2 

Now we can integrate over a, b' , w, and v, and arrive at 

expj^^Ny^ - N\Co\y2] 

(Cll) 
is obtained by substituting s = 



P(Co,Ci)cx / dy 



The final result 
lNy\ 



APPENDIX D: DISTRIBUTION OF Ai FOR /3 = 1 

In the large-A^ limit the joint distribution function 
P{A q,Ai) — P{Aq)P{Ai) factorizes as explained in Sec- 
tion |ni^ The distribution of Aq is given in Eq. (^. It 
remains to calculate the distribution of Ai — "Y^^TiUiVi. 
The Mi's and u^'s are independent Gaussian random num- 
bers. Averaging over them, we obtain the characteristic 
function 



x{p) 



(exp[iRe(pAi)]) 



\pn\ 



4iV2 
detiW'^Wf 



,det[(W^tM^)2 + |p|2/47v] 



(Dl) 



where p is a complex number. The Laguerre ensemble is 
again represented as the eigenvalues of the matrix prod- 
uct W^W, where W is the complex symmetric matrix 
with distribution (Al). Following the route of Appendix 
H we represent the determinant in the denominator by 
a Gaussian integral over a complex vector z and choose 
a basis in which W is of the form (|A8|). The charac- 
teristic function is then obtained as the following multi- 
dimensional integral. 



X(p) = /dxdydzd^ad^fody 



X |dety|4|a[6-r(r-i)ii]'-a;^|^exp - 



\zp\ 
4iV2 



X exp [~z^ {i\af + xY +x^\a + b*f + x^y^)] 

X exp [-jN{\a\^ + \b\^ + 2.t^ + V + tr Y^Y)] . 

(D2) 

Let us briefly describe in which order the integrations 
are performed most conveniently. Fourier transformation 
with respect to p converts the characteristic function back 
into the distribution function P{Ai). This step gives rise 
to a factor z~^ exp{—\Ai\'^N^z~^). We can also integrate 
over y, which results in a factor exp[— 27A^/(xz)^]. We 
introduce new variables by the substitutions b = b — a*, 
X = v/z, a ~ a! jz. After these transformations one 
succeeds to integrate over 6', z, and a' . The remaining 
integral over v = |v| is of the form 

P(Ai)cx|yli|-5 /df-y-^e-^/^ 



TT[8\Aif{2 + Av + v^) + 3^2(16 + 16t) + 3v^)] 
2\Ai\v 



[\Ai\''v^{288 + 30Av - 25v^) 



(|Al|2-K 1-2)4' 

4- |^i|%^(192 + 176i; - ITw^) + 8|Ai|*(6 + Av - w^) 
+ 3w*(16 + 16w + 3v^) + \Ai\^v^{l92 + 208w + Alv^)] 
- [16|Aip(2 -I- 4u + v^) + 6^2(16 + 16t; + Sw^)] 

V 



X arctan 



l^il 



(D3) 



The more compact form ( pq ) is the result of the replace- 
ment V ~ 2/s, followed by a number of partial integra- 
tions. 
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